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I Introduction 



Semileptonic weak hadron decays provide useful information on form factors of the 
weak current. The lepton pair can then be replaced with a hadron, permitting the 
calculation of nonleptonic decays. Although this hadron can re-interact with the 
rest of the system, the effects of this re-interaction sometimes can be neglected or 
evaluated. In such cases one is employing the factorization hypothesis. An early 
version of this hypothesis |l| was recently justified for certain decays of hadrons 
containing heavy quarks 0. 

In the present paper we update and test some factorization predictions first made a 
number of years ago We compare values of the Cabibbo-Kobayashi-Maskawa ma- 
trix element \V c b\ and form factor shapes obtained from (1) the differential distribution 
dT(B° -> D^+rpfi/dq 2 |, 3 g ft g g and (2) the color-favored two-body nonlep- 
tonic decays 5° — > D^ + (tt^, p~ , ajf). We find that consistency between nonleptonic 
and semileptonic determinations is at least as good as that among the semileptonic 
determinations themselves. 

Using a combined fit to semileptonic and nonleptonic decays and the measured 
value of the D s meson decay constant /d s , we then predict the rates for 5° — > 
£)(*)+£)(*)- and find that they are in accord with experiment. We thus find that fac- 
torization holds not only in color-favored cases in which the current produces a light 
meson, where it has been justified ||, but also when the current produces a heavy me- 
son, where no such justification has been presented. The importance of such processes 
has recently been stressed by Lipkin l lQfl . We also find that new experimental ratios 



of helicity amplitudes for color-favored processes agree with predictions. We shall 
ignore small non-factorizable contributions to color-favored B decays as discussed, 
for example, in Ref. [IT]. 



In Section II we review factorization predictions for the decays B° — > W*~ , 
where the virtual W*~ produces either a lepton pair l~Di or a hadron n~, p~ , aj~, 
D~, or D*~. These processes are purely color- favored. We do not consider the corre- 
sponding B~ decays, for which the nonleptonic processes receive both color-favored 
and color-suppressed contributions. We then (Section III) discuss the differential dis- 
tributions dT{B — > 1~ v{) / dq 2 and the information they can provide regarding 
the values of \V c b\ and the form factor slope at the normalization point. Results of fits 
to B° two-body decays to charmed final states are presented and compared with those 
from semileptonic decays in Section IV. We discuss the predictions of the factorization 
hypothesis for decays in which the weak current produces a in Section V and 
for ratios of helicity amplitudes in Section VI. Section VII concludes. An Appendix 
summarizes parameters of error ellipses used in combining data. 

II Notation and predictions 

We review notation which is described in more detail in Ref. ||. We consider processes 
in which a semileptonic 5° decay of the form shown in Fig. [l|(a) can be related to 
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Figure 1: Feynman diagrams for B decay illustrating the factorization hypothesis. 
A semileptonic decay (a) is related to a hadronic decay (b) by the replacement of the 
lepton pair of 4- momentum q with a qq' pair of effective mass q 2 . 



the corresponding decay with the lepton pair replaced by a quark pair, illustrated 
in Fig. |l](b). The matrix element for production of a pseudoscalar meson P(q) of 
4-momentum q from the vacuum by the axial vector current is 

(P(q)\A,\0) = tf P q, , (1) 

while that for production of a vector meson by the vector current is 

(V(q)\V fl \0) = e;M v f v , (2) 
and that for production of an axial vector meson by the axial vector current is 

(A(q)\A,\0) = e;M A f A , (3) 

The form factors for the ~B°(v) — > D^*\v') transitions are described in the heavy- 
quark limit by one universal function of the Lorentz-invariant variable w = v ■ v', 
where v and v' are invariant four- velocities: v = p^o/m B , v' = p D (*) /m D {*). We take 
c = 1 and note that q = p — p'. Another expression for w is then 

Pb-PdM m 2 B + m 2 M -q 2 
w = = . (4) 

m B m D {*) 2m B m D{ ,) 

[A variable (v — v') 2 = 2(1 — w) was called w in Ref. 0.] If e denotes the polarization 
vector of the final D*, we may write |l| 

(D(v')\V„\B(v)) = ^b^F v {w){v + v')^ , (5) 

(D*(v\ e)\A^\B(v)) = y/nw^F A (w)[(?(l + v ■ v') - e* • v v'^} , (6) 

(D\v\e)\V^\B{v)) = -i^Kb^F v {w)e^ a ^v a v' P . (7) 

We take m B = 5.28 GeV and = m D *+ = 2.01 GeV or m D + = 1.87 GeV depending 
on the final-state charmed meson. The maximum momentum transfer occurs when 
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the recoiling is at rest in the B rest frame, so g^ ax = (mg- m D (,)) 2 and hence 



w > 1. Another useful relation is 

,2 



= JL = i_ 2 um/c« + C w , c 



(*) _ m DM 



m 



B 



m. 



The differential decay width as a function of w for B° — > D^ + l v\ can then be 
written as 

^ = ^|K b | 2 m|m^ ) v / ^l/ W H|F y , A H| 2 , (9) 
where for B — ► -D/z// 



for B -> L>* Zi7 z 
and for £? — >■ 



/W = (™ 2 -i)(i±VC) 2 

/£H = 4u;y(u; + 1) , 



fin = a - vc*)> + ir • 



We shall take the form factors -FVa(w) to be parametrized as in Ref. The 



form factor i*V(iy) governing the process "i/j can be expressed as 



1 - 8p| y z + (51p 2 y - W)z 2 - (252p 2 Fv - 84)z J 



(10) 



where z = (y/w + 1 — \/2)/(\fw + 1 + v2) (the corresponding variable for elastic 
B ^ B transitions, a natural one for discussing analyticity in dispersion relations, 
was introduced in ||13||), while the form factor Fa(w) governing B — > is related 

to the axial- vector form factor A\(w) by 



1 + 



l-2wVC + C 

*\2 



1 + 



4w l-2wy/C + C 

W + 1 (1 - v^) 2 

u; — 1 



w + 1 



Ri(w) 



+ 



1 + 



F A (w) 
w — 1 



i-i? 2 H) 



can similarly be parametrized as 

= A x (l) x [l - 8pi^ + (53pi l - 15)^ 2 - (231^ - 91)z 3 



(12) 



These forms are motivated by dispersion relations |T2], [13], [14]]. R\(w) and R2(w) are 
given by 

i?iH = i?i(l)-0.12(w-l) + 0.05(«;-l) 2 , 

# 2 (w) = i? 2 (i) +Q.U(w- 1) -0.06(w- l) 2 . (13) 

In this paper we use the CLEO experimental results for Ri(l) and #2(1) M'- 



iJi(l) = 1.18 ±0.30 ±0.12, 
i? 2 (l) = 0.71 ±0.22 ±0.07. 



(14) 
(15) 



4 



As we know, p 2 F and p 2 Ai are the slope parameters for the form factors Fy(w) and 
Ai{w) at zero recoil, respectively. The difference between p Ai and p 2 F ■ , the slope 
parameter for F A {w) at w = 1, is predicted to be p F — p Al = 0.21 [T2|| . Recall that 
p 2 Fy and p 2 p A are actually the same in the single pole model. To make a connection 
between F v<A (w) and the single pole form factor |Q: 

Fv,aH = ^V,a(1)/[1 - 2(1 - ^)/^ 2 (yi A)], (16) 

we assume p|y = p 2 Ai —0.21. From now on we will simplify p\ as p 2 . This parameter 
describes the slope of the Isgur-Wise form factor at the zero-recoil point: p 2 = 
[dF Al (w)/dw]\ w=1 (see, e.g., PI). 

At w = 1, the vector and axial vector form factors are expected to behave as 
F V {1) = Vv (l + 8 1/mb ), F A (1) = VA (1 + 5 1/ml ). Here Vv = 1.022 ± 0.004 and 
Tj A = 0.960 ± 0.007 are QCD corrections |T7|. The terms <5i/ mb and S 1 / m 2 are non- 
perturbative in origin, and correspond physically to the excitation of states other 
than D and D*. Lacking a reliable method for estimating the term 5i/ mb in Fy(l) 



T8| , we set it equal to zero. The absence of a 5i/ mb term in F A (1) is the subject 
of Luke's theorem [l|. We take <5 1/m a = -0.05 ± 0.035 |2|, H, with the product 



Fa(1) = 0.913 ±0.042 as used in Ref. ||. Eq. (g) can then be integrated with respect 
to w to yield predicted decay rates as functions of the two parameters p 2 and \V c b\. 

The decay widths of some nonleptonic modes may be obtained under the assump- 
tion of factorization. For simplicity we assume all B° —>■ D + M~ transitions to involve 
Fv{wm) and all B — > D* + M~ transitions to involve F a (w1j), where 

w M = (m 2 B + m 2 D -m 2 M )/(2m B m D ) , w* M = (m 2 B +m 2 D * -m 2 M )/(2m B m D *) . (17) 



We then find 



G 2 



T(B°^D+n-) = ^\VMVu^mlfl\a x \ 2 \F v {w^)\ 2 {l- Jtf 



32tt' 



r(B° -> D* + n-) = ^I^HK.pm^laiHFAMI^l + A/C 7 )' 



32tt 1 



r(B°-Z>V) = #l^riK,rm|/ p 2 |a 1 | 2 |F y K)| 2 (l + ^ 



32tt 



xA^(i, C , C )MM^) (19) 



(20) 



xA i / 2 (i,c,g^-^ 



32tt' 



xA^CCp) ^^ (21) 



r(5°-D\) = ^|K b | 2 |K d | 2 m|/ 2 i |a 1 | 2 |F F («; ai )| 2 (l + V / C) i 



5 



xA^(l )C)Coi )^^ 



T(W -> D* + a~ 



where Cm = m M/ m B while 



G 2 



32tt 



:A 1/2 (l,C,C ai ; 



4\/C 



(22) 



(23) 



N(C, Cm) = (1 - VC*r[(l + VC*) 2 - Cm] + 4Cm(1 + C* - Cm) 
A(a, 6, c) = a 2 + 5 2 + c 2 - 2ab - 2ac - 26c . 



(24) 
(25) 



The QCD correction |ai| is taken to be 1.05 for all processes; this is a sufficiently good 
approximation to the actual situation, in which values differ by less than a percent 
from process to process 0. In the limit of small m^, the results ( [Op and ( |T9| ) are 
special cases of the simple Bjorken relation PJ 



r(5°^ J DW + 7 r-)=6 7 r 2 / 2 |K,r|a 1 | 



dq 2 



(26) 



III Semileptonic decays 

The CLEO Collaboration |§ at the Cornell Electron Storage Ring (CESR) and the 
ALEPH, DELPHI, and OPAL Collaborations §, |, [| at LEP gf have measured the 
spectra in lepton-pair squared effective mass q 2 (equivalently, in the Isgur-Wise vari- 
able w) for the decay B — > D* + l~i>i. The spectra may then be fitted for ^(1)1^1 
and p 2 . There is a strong correlation between the two parameters. The results 
are shown in Fig. |2| and Table [I], where we have taken Fa{1) = 0.913 ± 0.042 as 
in Ref. |3j]. Our fitted parameters for the CLEO data differ slightly from those 
presented by the collaboration itself, since we wished to generate an error ellipse 
and therefore fitted the spectral points directly without taking account of point-to- 
point correlations. For comparison, CLEO quotes \V c b\ = 0.0464 ± 0.0020(stat.) ± 
0.0021(syst.) ± 0.0021(theor.) and p 2 = 1.67 ± 0.11 ± 0.22. The combined fit implies 
\Vcb\ = 0.0399 ± 0.0023 and p 2 = 1.27 ± 0.26. 

The CLEO Collaboration has also measured the spectrum for the decay S° — > 
D + l~i?i P|. A fit to this spectrum with the form factor (10) and with Fy(l) = 1.022 
yields \V cb \ = 0.04S9±g;gg^ and p 2 = 1.33±g;|^. 



IV Nonleptonic two-body decays and combined fit 



We fit rates for nonleptonic two-body decays [Eqs. (|I8|-|2"3"|)1 to experimental averages, 
allowing (as in the fit to semileptonic spectra) for variation of \V c b\ and p 2 . We use 
the averages of Ref. [2~2"| except for modes for which new values have been presented 
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Figure 2: Error ellipses corresponding to A% 2 = 1 for fits to B — ► D* + l~i>i spectra. 
Dot-dashed lines correspond to fits to individual experiments, not including theoret- 
ical error in Fa(1)- Fit to combined semileptonic data without (with) error in Fa(1) 
is shown by the dashed (solid) ellipse. The errors we quote on each variable in Table 
| correspond to ±lcr extremes in which the other variable is also permitted to vary. 
The plotted points (+ for individual experiments, x for combined data) correspond 
to x 2 minima. 



7 



Table I: Values of | V c b\ and p 2 obtained from fits to individual spectra in B 
decays. 



Experiment 



cb\ 



CLEO i 
ALEPH § 

DELPHI § 
OPAL § 

Combined fa 



0.0461 ± 0.0036 
0.0361 ±0.0029 
0.0378 ± 0.0031 
0.0414 ±0.0023 
0.0399 ± 0.0023 



1.58 ±0.27 
0.74 ±0.50 
1.22 ±0.41 
1.44 ±0.32 
1.27 ±0.26 



(a) Errors include common theoretical error on Fa(1)- 



Table II: Branching ratios for B decays averaged for our fits, in units of 10" 



Mode 



Value 



Value 



Average 



D* + ti- 
D* + p- 
D* + D- 

d* + d:~ 



2.76 ±0.21 
6.8 ±3.4 H 



2.9 ±0.3 ±0.3 ]2J 
11.2 ±1.1 ±2.5 



2.79 ±0.19 
9.5 ±2.1 
10.0 ±3.1 

18.2 ± 3.7 ± 2.4 ± 4.6 H 25.4 ± 3.8 ± 5.3 ± 6.4 |5| 20.5 ± 6.3 



11.0 ± 1.8 ± 1.0 ± 2.8 gj 7.1 ± 2.4 ± 2.5 ± 1.8 |5 



p3| , [24] , |25| ; these are summarized and averaged in Table ([]]. For the decays involving 
(to be discussed in the next Section) the last errors in the second and third 
columns of Table [TJ refer to a common systematic error in D s branching ratios, which 
are based on B(D+ — > (f>7i + ) = (3.6 ±0.9)%, and are combined accordingly. In our fits 
we use \V U d\ = 0.974, r B = 1.548 ps, and f n = 131 MeV. f p and f ai are determined 
to be 209 MeV and 229 MeV ||, respectively, from the branching ratios for r — > pv 
and r — > a\u. 

The fit to two-body nonleptonic decays alone gives rise to a different correlation 
between |V^,| and p 2 than that to the B — > D^ + l~i>i spectra, since the decay rates 
are dominated by low q 2 and hence high w. Contours of A% 2 = 1 (la) for nonleptonic 
decays are shown along with the A% 2 = 1 contours for S° — > D^ + l^vi spectra in 
Fig. |3|. Also shown are the Ax 2 = 1 contours for the combined fit without and with 
common theoretical errors. We find \V c b\ = 0.0415±0.0022 and p 2 = 1.52±0.11. The 
results are summarized in Table III . The error on \ V c b\ is dominated by the theoretical 
uncertainty on the form factors at w = 1, which we take to have the same fractional 
value (0.042/0.913) for vector and axial form factors. 

The fitted branching ratios are compared with experimental data in Table [TV]. We 
also show the predictions of a recent investigation based on a more detailed application 
of the factorization hypothesis 0. The quality of the fit is acceptable except for a 
slight excess in the predicted branching ratio for S° — > D* + n~ . 

It is interesting to compare the form factors based on Eqs. (|T0D and ( |i~T|) with 
the simple pole model (|16|) [[|, where m B w (v,A) ^ as ^he interpretation of the mass 
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Combined Semileptonic and Nonleptonic Decays 
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Figure 3: Right-hand pair of dash-dotted curves: Contours of A% 2 = 1 (lcr) for fit 
to nonleptonic two-body decays B° — > D^ + (tt~, p~, aj~). Upper dash-dotted ellipse: 
Contour of A% 2 = 1 (lcr) for fit to B° — > D + l~i?i spectrum. Lower dash-dotted ellipse: 
Combined fit to semileptonic B — > D* + l~i?i decays. All these fits are performed 
assuming F A (1) = 0.913 and F v (l) = 1.022. Contour of Ax 2 = 1 (la) for fit to 
combined semileptonic and nonleptonic data without (with) common theoretical error 
in form factor normalization is shown by the dashed (solid) ellipse, x 2 minima are 
indicated by + for nonleptonic and B° — > D^ + l~V[ decays and by x for combination 
of all data. 



9 



Table III: Values of \V c b\ and p 2 obtained from fits to nonleptonic two-body B 
£)(*)+ (p- ? p- ; a -) decays and B° — > D^ + l~i?i decays. 





Decays 


\Vcb\ 


P 2 


B° 


Nonleptonic 


0.0450 (a) 


1.69 (a) 


— > D* + l~i?i spectra 


0.0399 ± 0.0023 


1.27 ±0.26 


B°- 


-> D + l~vi spectrum 


0.0459^;^! (a) 


1.33jl j5 ( a ) 




Combined (b) 


0.0415 ±0.0022 


1.52 ±0.11 



(a) Large correlated errors; see Fig. [|. 

(b) Errors include common theoretical error of 5Fa(1)/Fa(1) = 4.6%. 

Table IV: Branching ratios in units of 10 -3 : comparison between data and predictions. 



Decay mode Data Ref. || (a) Our fit x 2 contribution 



T? 


-> D + 7T- 


3.0 ±0.4 


3.27 


3.19 


0.22 


B° 


-> D*+TT- 


2.79 ±0.19 


3.05 


3.10 


2.70 


B° 


-> D + p- 


7.9 ±1.4 


7.64 


7.92 


0.00 


S° 


-> .D* + p- 


9.5 ±2.1 


7.59 


8.78 


0.12 


B° 


-> D+ar 


6.0 ±3.3 


7.76 


9.10 


0.88 


B° 




13.0 ±2.7 


8.53 


12.2 


0.09 


(a) 


For preferred values of form factors and = 1.05. 





of a pole in the weak b — > c (V,A) current. The axial form factor for p\ = 1.52 is 
compared with a pole model with wqa = 117 in Fig. |j. Also shown are CLEO [|J 
and DELPHI |J data points. The pole-model form factor is almost indistinguishable 



from that [0] motivated by dispersion relations. The value Wqa = 1.12 ±0.17, found 
in Ref. 0, is consistent with the present determination w a = 1.17 ± 0.08. 

The vector form factor (|T0|) is characterized by a slope parameter p 2 Fv = p 2 Ai —0.21 
i~2| and hence p 2 F = 1.31. It is compared to a pole form factor with wov = 1.14 and 
to the CLEO data in Fig. |5|. Thus, a nearly universal pole position characterizes 
the vector and axial form factors, as in Ref. 0. The CLEO data lie slightly above 
the predicted form factor but have the same w dependence, as one can also see in 
Fig. 3. It should be recalled that, in contrast to the case of the axial-vector form 
factor, theoretical estimates of the 0(1/ m&) correction to the vector form factor are 



lacking [18]. The normalization of the CLEO data may reflect our ignorance of this 



correction. 
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Figure 4: Form factors (|ll|) (solid curve) for p\ = 1.52 and pole-model (|T6|) (dashed 
curve) for wqa = 1.17. Data points are from CLEO (crosses) and DELPHI (dia- 
monds) . 
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Figure 5: Form factors fllCf) (solid curve) for p 2 Fy = 1.31 and pole-model (|I6| ) (dashed 
curve) for wqv = 1-14 compared with CLEO data (plotted points). 
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Table V: Comparison of predictions for branching ratios (in units of 10 3 ) involving 
production by the weak current in B° decays. 



Decay mode 


Data 


Prediction 


D+D- 


8.0 ±3.0 


14.9 ±4.1 


D* + D- 


10.0 ±3.1 


8.6 ±2.4 


D+D*- 


10.0 ±5.0 


10.0 ±2.8 


D* + D*- 


20.5 ±6.3 


24.0 ±6.7 



V Dy' production by the weak current 

When the qq' meson in Fig. Mb) is a D~ or D*~ , Eqs. ([T8|-p3D can be used for the 
respective predictions for r(2T -> D + D~), T(B° -> D* + D;), T(B° -> D + D*~), and 
T(B° -> D*+D*-) by replacing / w -> / Ds , / p -> -> Cc s = ™D s / m B> C P -> 

Co* = m 2 D ,/m 2 B , and other corresponding substitutions of kinematic variables u> and 
iV. We shall assume fz>* = fo s = 270±16±34 MeV based on an experimental average 
of rates for D s — > \w and D s — > rv [EB]; the latter error is the common systematic 
error associated with a 25% uncertainty in the branching ratio for Df — > 07r + . (For 
comparison, the value of 259 ± 74 MeV was found in Ref. || by utilizing the observed 
branching ratios of processes in which the weak current produced a D s or D*.) 

We then predict the branching ratios for production by the weak current 
shown in Table |V|. Experimental values are from Ref. [EjJ (D + D^~) or Table [TI| 
(D* + D^~). The predictions have an overall 25% uncertainty associated with all D s 
branching ratios. They are as well obeyed as those for the light mesons. 

An additional prediction involving heavy meson production by the weak current 
§ is that B(B° -> D*+D*-)/B(B° -> D*+D~) = 0.13(f D /f Da ) 2 ~ 0.09, where f D 
and fo a are the decay constants for the nonstrange and strange D mesons. The 

hO.04 
-0.03- 



experimental value for this ratio WJ\ is 0.067°'" 



VI Ratios of helicity amplitudes 

The decays of spinless particles to two vector mesons are describable by ampli- 
tudes A (longitudinal polarization), A\\ (linear parallel polarization) and A± (linear 
perpendicular polarization), normalized such that \A \ 2 ± \A\\\ 2 ± \Aj_\ 2 = 1. Factor- 



ization predicts (\A \ 2 , \A»\ 2 , \A ± \ 2 ) = (88, 10,2)% for B° -> D* + D*- and (55,39,6)% 







fori? — > D* + p . Experimental values are only quoted for | A \ 2 : (87.8 ± 3.4 ± 3.0)% 
for B° -> D* + D*- [E7| and (50.6 ± 13.9 ± 3.6)% for B° -> L>*+p- |29|. These agree 
with the predictions, as does the intermediate case of p'(1418) production |30 . 
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VII Conclusions 



New data on B meson decays have improved the precision of tests of some early 
factorization predictions |§, and yield a value \V c b\ = 0.0415±0.0022 when CLEO and 
LEP data on B — > D^ + l~ui spectra are combined with two-body nonleptonic decays 
B — > D^ + (ir~, p~, cli). The slope of the universal Isgur-Wise form factor at the 
normalization point w — 1 is found to be described by the parameter p 2 = 1.52±0.11. 
These values are only slightly different from those based on B — > D^ + l~i>i spectra 
alone, indicating that factorization for color-favored decays and universality of B — > 
£)(*) form factors are reasonable approximations. Consistency between nonleptonic 
and semileptonic determinations is at least as good as that among the semileptonic 
processes themselves. Our neglect of 0(1/ m^) corrections to the vector form factor 
may underestimate the rate for B° — > D + l~i>i slightly. 

Satisfactory rates for processes involving production by the weak current are 
obtained when the world average of direct measurements for fu g is used, and when it 
is assumed that fa* = fi> s . Ratios of helicity amplitudes for color-favored processes 
are also found to be in accord with predictions. 
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Appendix: Parametrization of error ellipses 

In order to combine results of fits to semileptonic decays in the absence of the ALEPH 
and DELPHI raw spectra, we have parametrized their fits in terms of error ellipses 
corresponding to A% 2 = 1, and generated corresponding ellipses for our own fits to 
CLEO and DELPHI data. The equations describing these ellipses are given below. 
Also shown are the contributions to A% 2 for each set of B — > D* + l~j?i data when 
x = p 2 and y = \V c b\ are taken to equal their values x = 1.52 and y = 0.0415 in the 
global fit. The corresponding A% 2 value is 4.02 for the sum of the nonleptonic modes 
listed in Table IV, which, when combined with the A% 2 values for the semileptonic 
spectra, leads to a total of A% 2 = 23.2 for the fit to five semileptonic spectra and six 
nonleptonic decay rates. The largest source of this Ax 2 (10.3) is the higher overall 
scale of the B° — > D + l^ui spectrum measured by CLEO, with some contribution also 
from the disparity between the CLEO and ALEPH fits. In comparison with these 
disagreements among purely semileptonic processes, the fits to nonleptonic decays do 
not fare badly at all. 
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CLEO: 



41.8814(x-x c ) 2 -5174.3(x-x c )(y-y c )+238793(?/-?/ c ) 2 = 1; 

x c = 1.58; y c = 0.0461; A X 2 = 3.78; 



ALEPH: 



9.6816(a; - x a f - 2554.0(x - x a ) (y - y a ) + 287560(y - y a f = 1; 

x a = 0.74; y a = 0.0361; A X 2 = 3.52; 



OPAL: 

19.4601(x-x o ) 2 - 3771.3(x-x o )(y-?/ o ) + 379980(y-y o ) 2 = 1; 

x Q = \M; y Q = 0.0414; A^ 2 = 0.10; 

DELPHI: 

14.1777(x- a;,) 2 - 2856.4(x - a; d )(y - y d ) + 246120(y - y,) 2 = 1; 

x d = 1.22; y d = 0.0378; A X 2 = 1.47; 

CLEO+ALEPH+OPAL+DELPHI (without common theoretical normalization 
error; divide all coefficients by 2.850 for ellipse with this error): 

85.2008(x-x) 2 - 14356(x-x)(y-y) + 1152453(y -y) 2 = 1; 
x = 1.27; y = 0.0399; 

CLEO+ALEPH+OPAL+DELPHI+DECAYS (without common theoretical nor- 
malization error; divide all coefficients by 3.514 for ellipse with this error): 

1059.7(x-x) 2 - 86781(x-x)(y-y) + 252U00(y-y) 2 = 1; 
£ = 1.52; y = 0.0415; 
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